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The quantum kicked rotator is a periodically driven dynamical system with a metal-insulator 
transition. We extend the model so that it includes phase transitions between a metal and a 
topological insulator, in the universality class ol the quantum spin Hall effect. We calculate the Z2 
topological invariant using a scattering formulation that remains valid in the presence of disorder. 
The scaling laws at the phase transition can be studied efficiently by replacing one of the two spatial 
dimensions with a second incommensurate driving frequency. We find that the critical exponent 
does not depend on the topological invariant, in agreement with earlier independent results from 
the network model of the quantum spin Hall effect. 
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I. INTRODUCTION 

The spin of an electron moving in an electric field expe- 
riences a torque, which can be understood as arising from 
the magnetic field in its rest frame. In a two-dimensional 
electron gas this velocity-dependent magnetic field pro- 
duces the quantum spin Hall effect, 1,2 reminiscent of the 
quantum Hall effect but without time-reversal symme- 
try breaking. 3 The difference manifests itself in the inte- 
ger quantized values Q of the dimensionless conductance. 
While there is no restriction on Q 6 Z in the quantum 
Hall effect (QHE), only the two values Q = 0, 1 e Z 2 
appear in the quantum spin Hall effect (QSHE). In both 
effects the current is carried by edge states, separated by 
an insulating bulk. The insulator in the QSHE is called a 
topological insulator 4 ' 5 if the topological quantum num- 
ber Q = 1 and a trivial insulator if Q = 0. 

Bulk states delocalize when the conductance switches 
between quantized values. A distinguishing feature of the 
QSHE is that the delocalized states can support metallic 
conduction (conductance e 2 /h), while in the QHE the 
conductance remains < e 2 /h. The metallic conduction 
appears in extended regions of phase space, separated by 
a quantum phase transition (Anderson transition) from 
the regions with a quantized conductance and an insulat- 
ing bulk. For a trivial insulator this is the familiar metal- 
insulator transition in a two-dimensional (2D) electron 
gas with spin-orbit scattering. 6 ' 

A fundamental question raised by the discovery of the 
QSHE was whether the phase transition would depend 
on the topological quantum number. Specifically, is the 
critical exponent vq of the diverging localization length 
different if the phase boundary separates a metal from 
a topological insulator, rather than a trivial insulator? 
A numerical simulation 8 of the Kane-Mele model 1 of 
the QSHE gave an affirmative answer, finding a value 
v\ « 1.6 substantially below the established result 9 ' 10 for 
the metal-trivial-insulator transition. In contrast, stud- 



ies of the network model in the QSHE universality class 
gave v\ = vq sa 2.7 within numerical accuracy, 11,12 con- 
sistent with analytical considerations from the nonlinear 
sigma model 13 why the critical exponent should be Q- 
independent. 

In this work we study the metal-insulator transition 
in the QSHE by means of an altogether different, stro- 
boscopies model — the periodically driven system known 
as the quantum kicked rotator. 14-16 A key feature of this 
dynamical system is that spatial dimensions can be ex- 
changed for incommensurate driving frequencies, 17,18 al- 
lowing for the study of metal-insulator transitions in one 
spatial dimension. 19,20 This proved very effective in the 
QHE, 22 and also made it possible to experimentally study 
the 3D Anderson transition in a ID optical lattice. 21 In 
the present paper we apply the same strategy to study 
the 2D QSHE in a ID system. 

In the next section we show how the quantum kicked 
rotator can be extended to include the topological Q = 1 
phase of the QSHE. We first construct this Z2 quantum 
kicked rotator in 2D and then carry out the mapping to 
ID. We calculate the phase diagram in Sec. Ill, using a 
scattering formula for the topological quantum number 
valid for disordered systems. 23 ' 24 In Sec. IV we determine 
the scaling law at the metal-insulator phase transition 
and compare the critical exponents vq for Q = and 
= 1. We conclude in Sec. V. 



II. CONSTRUCTION OF THE Z 2 QUANTUM 
KICKED ROTATOR 



A. Stationary model without disorder 

We start from a translationally invariant 2D system 
and will add disorder later. The minimal model Hamil- 
tonian Ho(p) of the QSHE has four bands at each mo- 
mentum p = (pi,P2), distinguished by indices a (up and 
down spins) and r (s and p orbitals). The Pauli matrices 
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<Ji and Tj (i = 0,x,y,z) act on the spin and orbital de- 
grees of freedom, respectively. Time-reversal symmetry 
is essential, 

a y H*(-p)a y = H{p). (2.1) 

Inversion symmetry is not essential (and will be broken 
anyway once we add disorder) , but is assumed for conve- 
nience, 

r z H(-p)r z = H(p). (2.2) 

The generic Hamiltonian that satisfies the symmetries 
(2.1) and (2.2) has the form 25 

5 

H(p) = E (p) + Y, Up)T a , (2.3) 

a=l 

r = (t x (J z , Ty(TQ, T Z (T , r x a x , T x (Jy). (2.4) 

The real functions Eq, f$ are even under inversion of p, 
while the functions fi, /a, Ja, fb are odd. Because the T- 
matrices anticommute, {r a ,r^} = 28 a p, the band struc- 
ture is given by 

e±(p) = Eo(p) ± y/EJZfcj. (2.5) 

Each band is twofold degenerate. 

The band gap can close upon variation of a single con- 
trol parameter at high-symmetry points A a in the first 
Brillouin zone, which satisfy A a = — A a modulo a re- 
ciprocal lattice vector. At these time-reversally invariant 
momenta the Bloch wave function w_(A a ) of the lower 
band has a definite parity ir a = ±1 under inversion, 
r 2 u_(A a ) = 7r a u_(A a ). The Z 2 topological quantum 
number Q follows from the Fu-Kanc formula, 26 

(-i)«=n*«- (2.6) 

a 

A gap closing and reopening can switch the parity of 
the lower and upper bands, inducing a change in Q (a 
topological phase transition). 

The specific choice for the functions f a (p) which we 
will study in the following is based on experience with 
the stroboscopic model of the QHE. 22 We take 

£0 = 0, f a = T(u)u a , T(u) = (2/w)arctanw, (2.7) 

where the vector u (of length u = \u\) has components 

u {p) = (-K" sin Pii Ksinp2, /3K(fj, — cospi — C0SP2), 

7A^cospi sinp2, jKcosp2 sinpi). (2-8) 

For 7 = this is the Bernevig-Hughes-Zhang model of 
the QSHE, 27 up to a function T(u) > which flattens 
the bands without closing the band gap (hence without 
affecting the topological quantum number) . Without the 
7-term the spin degree of freedom a is conserved and 
the QSHE Hamiltonian is identical to two copies of the 
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FIG. 1: (Color online) Band structure for the clean Hamil- 
tonian H(p), as calculated from Eqs. (2.5), (2.7), and (2.8), 
for parameter values K = 2, j3 = 0.8, 7 = 2, [i = —0.3. The 
Dirac cones emerging at the time-reversally invariant points 
in the Brillouin zone will touch when fi = 0. 

QHE Hamiltonian (with opposite magnetic fields, to re- 
store time-reversal symmetry). When 7 7^ spin-orbit 
coupling mixes the spin-up and spin-down blocks of the 
Hamiltonian. 

The band structure for one set of parameters is shown 
in Fig. 1. The time-reversally invariant momenta A a in 
the first Brillouin zone — ir < Pi,P2 < tt are at the four 
points (0,0), (7r,7r), (0,tt), (w, 0). (We have set both h 
and the lattice constant a equal to unity.) The topologi- 
cal quantum number (2.6) depends only on the parameter 
Mi 

f if > 2, 
Q - { 1 if \fx\ <2. (2 ' 9) 

For /x = 2 or — 2 the band gap closes at (0,0) or (tt, n), 
with a switch in parity and a change of Q. For fi = the 
gap closes at both points (0, w) and (tt, 0) — at constant 
Q since the two parity switches cancel. 



B. Time-dependent model with disorder 

The time-dependent model is based on the quantum 
kicked rotator 14 , which is a dynamical system designed 
to study the localization by disorder with great nucrical 
efficiency 15 ' 16,20 The time dependent Hamiltonian H(t) 
contains a disorder potential V(x) and a stroboscopic 
kinetic energy H(p), 

00 

H(t)=V(x)+H(p) J2 s (t~ n )- ( 2 - 10 ) 

n— — 00 

(The stroboscopic period r has been set equal to unity.) 
We take H{p) of the form (2.3) and will specify V(x) 
later. The disorder strength is set by the relative impor- 
tance of V(x) and H{p). The Floquet operator T de- 
scribes the time evolution of the wavefunction over one 
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period, 



iV{id p ) 



(2.11) 



Here i(9p is the operator x in the momentum represen- 
tation. The eigenvalues e~ l£ of the unitary operator T 
define the quasi-encrgics e <E (— 7T, 7t). 

We use the 27r-periodicity of i? (p) to label the eigen- 
states 'iqip) of J 7 by a Bloch vector q in the Brillouin 
zone — 7r < qi, <j2 < 7r. By construction, 



with Xg(p) a 27r-periodic eigenstate of 



e -iff(p) e -iV(ia p +g) 



(2.12) 



(2.13) 



This is the quantum kicked rotator with Z2 symmetry. 



C. Mapping from 2D to ID 

By adding an incommensurate driving frequency to 
the quantum kicked rotator in two spatial dimensions, 
it is possible to simulate the system in one single 
dimension. 17,18 For that purpose we take a potential of 
the form 



V(x) = Vi(xi) - m 2 , 



(2.14) 



with w/27r an irrational number G (0,1). We consider 
states which at t = are plane waves in the ^-direction, 
having a well-defined initial momentum p 2 = a. In the 
Floquet operator the term linear in 22 has the effect of 
shifting P2 to p 2 + w, so that the 2D operator T q can be 
replaced at the n-th time step by the ID operator 



x-(n) 
J q 



g— iff(pi, nu+a) - iVi(iO P1 +g) 



(2.15) 



There is considerable freedom in the choice of the poten- 
tial Vi in the remaining dimension. The simple quadratic 
form 



Vi(xi) = \{xi - .t )' 



(2.16) 



provides sufficient randomness if {w, A, 2-7r} form an in- 
commensurate triplet. 28 We take A = 1, 

The numerical simulation is performed by subsequent 
multiplication of the state Xg( x l 1 1) with the series of Flo- 



quet operators J 7 , 



O) 



0, 1,2,.. ., in the plane wave ba- 



sis e~ lmpi (eigenstates of x\). The integer m is restricted 
to the values 1,2, ... M, with M an even integer that sets 
the system size. As initial condition we choose 



(2.17) 



spatially localized at xq = M/2. The vector S is a nor- 
malized vector of rank 4 with random components denot- 
ing the spin and orbital degrees of freedom. 

The multiplication with the Floquet operators, repre- 
sented by M x M unitary matrices, can be done very 




FIG. 2: (Color online) Bottom panel: Time-dependent diffu- 
sion coefficient (3.2) as a function of /1 for K — 2, f3 — 0.8, 
7 = 2, shown for times t = 10 6 (red), t = 10 5 (green) and 
t = 10 4 (blue). The points of intersection of these curves lo- 
cate the metal-insulator transition. Top panel: Topological 
quantum number (3.3) for the same parameter values, used 
to distinguish the topologically trivial (Q — 0) and nontrivial 
(Q = 1) insulators. In the metallic regions Q is not quantized. 



efficiently by means of the Fast Fourier Transform algo- 
rithm. (See the analogous calculation in the QHE for a 
more detailed exposition. 22 ) The calculation is repeated 
for different values of a and q, to simulate a disorder av- 
erage. The disorder strength can be varied by varying 
K, with small K corresponding to strong disorder. 



III. PHASE DIAGRAM WITH DISORDER 

Before embarking on a calculation of the scaling law 
and critical exponents, we first locate the metal- insulator 
transitions and identify the topological phases. 

To find the metal-insulator transitions, we calculate 
the time dependent diffusion coefficient 



D(t) 



A 2 (t) 



t 



from the mean squared displacement 



A 2 (t) = ((*!(*) 



zo) 2 )- 



(3.1) 



(3.2) 



The brackets (• • • ) denote the expectation value in the 
state Xq( x iit) an d the overline indicates the ensemble 
average over a and q. We typically average over 10 
samples of size M = 4 ■ 10 3 . 

A representative series of scans of D versus fj, for dif- 
ferent values of t is shown in Fig. 2 (bottom panel). In 
Fig. 3 we show D versus t for different fi. In the insulat- 
ing phase D(t) oc 1/t decays with increasing time, while 
in the metallic phase D(t) oc ln£ grows with increasing 
time. 29 The metal-insulator transition at [1 = fx c is sig- 
naled by a t- independent D(t) (corresponding to a scale 
invariant diffusion coefficient). 
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FIG. 3: (Color online) Time-dependent diffusion coefficient 
as a function of t for several values of jj, near the metal- 
topological-insulator transition. (Same parameters as in Fig. 
2.) The slope vanishes at the transition point \i c ~ —0.368. 



In this way we can locate the phase boundaries, but we 
cannot yet distinguish topologically trivial and nontriv- 
ial insulators. For that purpose we need the topological 
quantum number. The formula (2.6) for the topological 
quantum number of the translationally invariant system 
does not apply for nonzero disorder potential. The scat- 
tering formulation 23,24 continues to apply and is what 
we will use. An alternative Hamiltonian formulation for 
disordered systems has been given by Prodan. 30 ' 31 

The topological invariant is computed using the for- 
malism described in Ref. 22. From the Floquet operator 
we can construct a reflection matrix r(s, <p) for a cylin- 
drical system enclosing a flux $ = <j>h/e. The topologi- 
cal invariant then follows from the following combination 
of determinants and Pfaffians, evaluated at e = and 

(j) = 0,7T, 

( pQ = %M Vdctr(0,0) 

1 ' Pf[(7„r(0,0)] ^detr(0,7r)' K ' } 

The results in Fig. 2 (top panel) show the disorder- 
averaged /^-dependence of Q for a system of size M x M — 
30 x 30. The value of Q is only quantized £ {0, 1} in the 
insulating regions. In the metallic regions Q averages to 
1/2 for a sufficiently large system, 23,31 which is not quite 
observed for our system sizes. 

The phase diagram obtained in this way is shown in 
Figs. 4 and 5. Without disorder, the topological invari- 
ant (2.9) gives a topological insulator for \fi\ < 2 and a 
trivial insulator for \fi\ > 2. With disorder a metallic 
phase appears between the insulating phases, provided 
that 7//3 0. For small 7 we find an additional triv- 
ial insulating phase around fi — 0, consistent with what 
was found in the quantum Hall system in the presence of 
disorder 22 . 

In some other models of the QSHE a reentrant behav- 
ior Q = 1— > 1 1— > with increasing disorder is observed 
in some regions of parameter space. 31-34 In our model a 
system which starts out topologically trivial in the clean 
limit stays trivial with disorder. 
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FIG. 4: (Color online) Phase diagram of the Z2 quantum 
kicked rotator in the fi-l/K plane for fixed ft = 0.8 and 7 = 2. 
The topological quantum number Q — 0, 1 distinguishes the 
topologically trivial and nontrivial insulating phases. The 
metallic phase between the insulating phases disappears in 
the clean limit 1/K — ► 0. The two circles indicate the metal- 
trivial-insulator and metal-topological-insulator transitions 
that were studied to compare the critical exponents. 



IV. SCALING LAW AND CRITICAL 
EXPONENT 

The premise of one-parameter scaling is expressed by 
the equation 35 

ln£>(t)=^(r a *), (4.1) 

where T is a universal scaling function. The localization 
length £ has a power law divergence oc \\x — ^ c \~ u at the 
metal-insulator transition, with critical exponent v. 

We follow the established method of finite-size (here: 
finite-time) scaling to extract v from the numerical 
data. 36 We rewrite the scaling law (4.1) as 

In D(t) = Fiut 1 / 2 "), (4.2) 
u= (n- /i c ) + ^ c fe (/i - fi c ) k , (4.3) 

k=2 

where F(z) is an analytic function of z = v£ L l 2v . By 
fitting the free parameters of the series expansion 

lnD(i) =hxD c + ^d k (ut 1 / 2 ") +S(t) (4.4) 

k=l 

to the data for D as a function of /1 and t, the critical 
exponent v is obtained. The extra term S(t) accounts 
for finite-time corrections to single-parameter scaling, of 
the form 

S(t)=t-yY, 9k ^-fx c )t 1 / 2 ") . (4.5) 

fe=0 

We have considered times up to f = 10 6 for system 
sizes up to M = 10 4 . The number of terms N u , N F , 
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FIG. 5: (Color online) Phase diagram for fixed disorder 
strength (K — 2), in the fi-y plane for f3 = 0.8 (top panel) 
and in the fi-/3 plane for 7 = 2 (bottom panel). The metal- 
lic phase is stabilized by increasing the spin-orbit coupling 
strength 7 or by decreasing /3. 





N u 


N F N s 


V 


X 2 /ndf 




2 


3 2 


2.67 ±0.09 


0.97 




2 


1 2 


2.69 ±0.06 


0.92 



TABLE I: Results of the finite-time scaling analysis described 
in the text. 



V. CONCLUSION 

In conclusion, we have presented a numerical method 
to study the metal-insulator transition in the quantum 
spin Hall effect (QSHE), based on the incorporation of 
Z2 topological symmetry into the quantum kicked rota- 
tor. We find that the critical exponent vq of the diverging 
localization length is the same whether the metal is ap- 
proached from the topologically trivial insulator (Q = 0) 
or from the topologically nontrivial insulator (Q = 1). 
Our results i>q = V\ 2.7 are in agreement with Refs. 
11,12, but not with Rcf. 8 (which found a much smaller 
v\ ~ 1.6 for the topologically nontrivial insulator). Since 
our stroboscopic model is fully independent of the net- 
work model used in Refs. 11,12, this is significant support 
for the insensitivity of v to the topological quantum num- 
ber Q. 

A special feature of the quantum kicked rotator is that 
it allows the study of the QSHE, as well as the QHE, 22 
in one spatial dimension, by exchanging a dimension for 
an incommensurate driving frequency. 17 ' 18 There is much 
interest in such one-dimensional models of topological 
phases, 37-40 because they might be more easily realized 
in optical lattices of cold atoms than the original two- 
dimensional models. Critical exponents were not stud- 
ied in these earlier investigations, which focused on the 
QSHE in clean systems without disorder. 

An interesting direction for future research is to study 
the edge states in the Z2 quantum kicked rotator, by 
replacing the periodic boundary condition used in this 
work by a zero-current boundary condition. While the 
localization length exponent v does not depend on the 
topological quantum number, the edge state structure 
does depend on Q, with a characteristic multifractality 
at the metal- insulator transition. 11 



and N s in the series expansions (4.3)-(4.5) are system- 
atically increased until the chi-square-value per degree 
of freedom (x 2 /ndf) is approximately unity, see Table I. 
The calculation is carried out at the two points indicated 
by circles in Fig. 4, one a metal-trivial-insulator transi- 
tion (giving vq = 2.67 ± 0.09) and the other a metal- 
topological-insulator transition (giving V\ = 2.69 ±0.06). 
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